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摘要 : 对 于 扎 商 甲骨 文 卜 辞 中 存在 的 唯一 可 考 的 一 次 日 食 的 记载 ， 作 者 利用 美国 
MEMA (NASA) 喷气 推进 实验 室 (Jet Propulsion Laboratory, JPL) 发 布 的 
DE422 行星 历 表 ， 解 算 了 记载 涉及 此 乙 已 日 食 的 主要 参数 ， 给 出 了 日 食 带 在 地 球 
上 的 覆盖 范围 。 进 一 步 确认 了 此 次 日 食 在 息 都 安阳 的 可 视 性 ， 验 证 了 有 息 卜 辞 对 此 
次 日 食 预 测 的 准确 度 。 
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我 国 是 世界 上 天 文学 发 达 最 早 的 国家 之 一 ， 早 在 三 千 多 年 前 就 有 可 考 的 文字 
星象 记载 。 这 些 古 代 天 象 记载 不 仅 对 于 历史 研究 和 古代 天 文学 研究 有 重要 价值 ， 
在 现代 天 文学 研究 中 也 起 着 重要 的 作用 由。 有 息 卜 辞 是 甲骨 文 的 主要 内 容 之 一 ,是 研 
究 商 周 社会 历史 的 重要 资料 。 谍 卜辞 曾经 记录 一 次 对 乙 已 日 食 完整 的 预测 ， 从 王 
子 日 开始 卜 问 ， 两 天 后 是 否 有 日 食 ? HAAN. SB. RR, YESH, 12 
次 不 厌 其 烦 地 卜 问 ， 直 到 第 54 日 的 乙 已 日 ， 史 官 记载 “日 有 食 ， 夕 告 于 上 甲 ， 九 
牛 ” 的 验 词 , 意思 是 : “日 食 观 测 到 了 ， 晚 上 为 此 告 祭 先王 上 甲 ， 献 祭 了 九 头 牛 。” 
这 也 是 息 商 甲骨 文 卜辞 中 唯一 可 考 的 一 次 日 食 预测 。 冯 时 进一步 综合 多 位 学 者 的 
研究 成 果 ， 并 根据 干支 的 记载 及 不 同时 期 的 息 王 所 处 的 年 代 ， 考 订 认为 此 卜辞 涉 
及 的 是 发 生 于 公元 前 1161 年 10 月 31 日 的 日 食 站 。 张 培 瑜 利 用 现代 天 文 年 历 依据 
的 纽 康 给 出 的 根 数 和 公式 ， 得 出 太阳 黄 经 ， 再 用 线性 内 插 分 析 了 公元 前 1500 年 到 
公元 2050 年 间 的 天 象 ， 并 给 出 了 中 国 古代 历史 名 城 的 见 食 情况 外。 刘 次 沅 和 马莉 
萍 借助 Kluepfel 的 SunTracker 软件 计算 日 食 参 数 ,给 出 公元 前 2300 年 到 公元 2100 
年 间 在 中 国 可 见 日 食 的 大 致 覆盖 范围 和 内。Espenak 和 Meeus 利用 VSOP87 和 
ELP2000/82 计算 日 月 位 置 ， 发 布 了 公元 前 2000 年 到 公元 3000 年 间 日 食 的 地 面 可 
视 情 况 59。 本 文 利用 美国 航空 航天 局 (NASA) 喷气 推进 实验 室 (Jet Propulsion 
Laboratory, JPL) 发 布 的 DE422 行星 历 表 计算 日 月 位 置 ， 研 究 公 元 前 1161 年 10 
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H 31 日 日 食 的 日 食 带 分 布 ， 重 点 考察 股 都 安阳 对 此 次 日 食 的 可 见 性 ， 为 深入 研究 
此 次 甲骨 文 日 食 的 记载 提供 天 文学 帮助 。 


1 天 文 行星 历 表 简 述 

天 文 行星 历 表 对 大 地 测量 、 天 体 测量 、 深 空 探测 、 飞 行 器 测控 及 其 他 相关 学 
科 领 域 均 都 具有 重要 的 实用 价值 。 根 据 编制 原理 ， 太 阳 系 行星 和 月 球 历 表 主 要 分 
为 数值 历 表 和 半 解 析 历 表 两 大 类 。 数 值 历 表 主 要 有 IPL 发 布 的 DE 系列 行星 历 表 、 
法 国 天 体力 学 与 历 算 研究 所 发 布 的 巴黎 天 文 台 行 星 数 值 积 分 系列 历 表 、 俄 罗斯 应 
用 天 文 研究 所 发 布 的 行星 和 月 球 系列 历 表 。 半 解析 历 表 主 要 为 法 国 天 体力 学 与 历 
算 研究 所 巴黎 月 球 历 表 (ELP2000/82) 、 行 星 轨 道 的 长 期 变化 序列 历 表 
(VSOP82/VSOP87/VSOP2000/VSOP2013) 、 外 行星 理论 系列 历 表 等 。 目 前 ，DE 
系列 行星 历 表 已 经 在 天 体 测量 、 深 空 导航 、 行 星际 探测 等 方面 得 到 较为 广泛 的 应 
H. 2009 年 9 月 JPL 发 布 的 DE422 行星 历 表 ， 数 据 跨度 自 公 元 前 3001 年 12 月 
7 日 至 公元 3000 年 1 月 30 日 ， 并 包含 章 动 和 天 平 动 信息 。 该 历 表 可 以 满足 研究 者 
对 古代 天 象 记载 研究 的 基本 需求 1。 


2 基于 DE422 行星 历 表 的 日 食 计 算 
2.1 日 月 的 视 位 置 

为 了 计算 日 食 的 有 关 参 数 ， 需 要 知道 在 地 心 处 观测 日 月 的 精确 位 置 。 进 入 20 
世纪 90 年 代 以 来 ，JPL 发 布 的 DE 系列 行星 历 表 都 采用 国际 天 文联 合 会 
(International Astronomical Union, IAU) 推荐 的 国际 天 球 参考 系统 international 


Celestial Reference System, ICRS) ， 对 应 的 国际 天 球 参 考 架 (International Celestial 
Reference Frame, ICRF) 是 通过 一 套 河 外 射电 源 的 位 置 来 实现 。 其 坐标 原点 是 太 
阳 系 质心 , 坐标 轴 的 指向 在 太空 中 是 固定 的 .DE422 行星 历 表 也 采用 ICRF 参考 架 。 
利用 该 历 表 可 以 直接 内 插 得 到 数据 跨度 内 行星 在 任 一 时 刻 的 三 维 直角 坐标 和 速度 
言 息 。 进 一 步 利用 美国 海军 天 文 台 提供 的 NOVAS 软件 ,可 以 得 出 日 月 的 视 赤 经 和 
IIRA o 
2.2 R hÈLEHA 

(1) 日 月 交会 时 刻 


本 文 考虑 公元 前 1161 年 10 H 31 日 发 生 的 日 食 ， 首 先 内 插 得 到 日 月 视 赤 经 相 
同 的 时 刻 为 地 球 动力 学 时 07:30:24 (Terrestrial Dynamical Time，TDT)， 这 也 是 日 


月 交会 时 刻 。 进 一 步 分 析 表 明 ， 此 次 日 食 为 日 环 食 。 表 1 给 出 发 生 初 亏 、 食 既 、 
食 甚 、 生 光 、 复 圆 等 食 相 的 TDT 时 刻 。 


#21 公元 前 1161 年 10 月 31 日 日 食 食 相 的 TDT 时 刻 
Table 1 TDT time of solar eclipse phase on October 31, 1161 BCE 


Partial eclipse Total eclipse Maximum Total eclipse Partial eclipse 
begins begins eclipse ends ends 
04:30:32 05:39:02 07:36:59 09:35:01 10:43:31 
(2) 日 食 带 分 布 


考虑 到 地 球 自转 的 长 期 变化 , 地 球 动力 学 时 (TDT) 和 世界 时 (Universal Time, 
UT) 之 间 存 在 一 个 时 间 改 正 数 AT" ，AT 数 值 随 时 间 发 生变 化 。 利 用 古代 的 天 
象 记载 情况 ， 可 以 推算 古代 的 AT 数值 。 反 之 ， 已 知 古代 某 天 象 年 代 的 AT 数值 ， 
可 以 推算 该 天 象 当时 的 可 观测 情况 。Stephenson 综合 古代 的 若干 天 象 记载 ， 设 定 月 
球 本 征 加 速度 为 -26.0 和 角 秒 /( 百 年 》>， 给 出 AT 的 分 段 拟 合 公式 中 1。 本文 所 采用 的 
DE422 行星 历 表 中 , 月 球 本 征 加 速度 为 -25.85” 人 (百年 > "4 加 入 该 差别 引入 的 改正 
后 ， 可 以 得 到 公元 前 1161 年 10 月 31 日 的 时 间 改 正 数 AT 为 27361 秒 。 根 据 日 月 
位 置 形成 的 日 食 几何 关系 上 ， 可 以 得 到 此 次 日 食 的 分 布 情况 〈 见 图 1) 。 
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. -1160.10.31 
Annular eclipse (JDo 1297671.5) 


Greatest Ecl. Time 


AH SS | 07:36:59 


Mag. =0.96 Wide =299(km) 
图 1 公元 前 1161 年 10 月 31 日 的 日 食 带 分 布 (方位 角 投影 ) 
Fig.1 The solar eclipse path on October 31, 1161 BCE under azimuthal projection 


图 1 中 Mag 为 食 分 ，Wide 为 中 心 食 的 日 食 带 宽度 。 图 2 以 圆柱 投影 方式 给 出 
此 次 日 食 的 日 食 带 分 布 。 
Annular eclipse -1160.10.31 (JDo 1297671.5) 


Greatest Ecl. t=07:36:59 AT= 27361 sec 
-150° -120° -90° -60° -30° 


Sa z x z ae ae ee ae 
Mag. =0.96 Wide= 299 (km) 


图 2 公元 前 1161 年 10 月 31 日 的 日 食 带 分 布 “圆柱 投影 ) 
Fig.2 The solar eclipse path on October 31, 1161 BCE under cylinder projection 
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为 了 直观 展示 古代 中 国 中 原 地 区 对 此 次 日 食 的 可 观测 情况 ， 对 图 2 局 部 放大 
得 到 图 3。 图 中 ， 息 都 安阳 (经度: 114°.38, AEE: 36.1) 以 黑色 五 角 星 标 出 。 


Annular eclipse -1160.10.31 (JDo 1297671.5) 
Greatest Ecl. t=07:36:59 AT= 27361 sec 


Wide= 299 (km) 


图 3 图 2 的 局 部 放大 


Fig.3 Partial enlargement of Fig. 2 


3 小 结 

本 文 利用 美国 喷气 推进 实验 室 发 布 的 DE422 行星 历 表 ， 计 算 了 有 扎 商 甲骨 文 小 
辞 中 唯一 可 考 的 公元 前 1161 年 10 月 31 日 乙 已 日 食 ， 给 出 日 食 带 分 布 ， 确 认 了 此 
次 日 食 在 笑 都 安阳 的 可 视 性 ， 此 次 日 食 在 中 国 中 东部 地 区 都 是 带 食 而 出 。 本 文 结 
果 表 明 ， 括 都 安阳 对 此 次 日 食 具 有 实际 观测 到 的 条 件 ， 从 侧面 验证 了 笑 下 辞 对 此 
次 日 食 预 测 的 准确 度 。 
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Research of Yisi Solar Eclipse Observation in the Oracle-bone 


Inscriptions of the Yin Dynasty Based on the JPL Planetary Ephemeris 

Ma Lihua, Han Yanben, Yin Zhigiang, Qiao Qiyuan 
(National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, 
Email: mlh @nao.cas.cn) 
Abstract: For the only solar eclipse recoded in the oracle-bone inscriptions of the Yin 
dynasty discovered so far, the authors use the DE422 planetary ephemeris released by 
the Jet Propulsion Laboratory (JPL) of the National Aeronautics and Space 
Administration (NASA) to calculate the main parameters involved in the Yisi solar 
eclipse, give the coverage of the solar eclipse path on the Earth. Furthermore, the 
visibility of the solar eclipse in Anyang, capital of the Yin dynasty is confirmed and the 
accuracy of the Oracle-bone inscriptions about the solar eclipse of the Yin dynasty is 
verified. 


Keywords: Oracle-bone inscriptions; JPL planetary ephemeris; solar eclipse path 


